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ABSTRACT 

We look for evidence of tidal stripping in elliptical galaxies through the analysis of homoge- 
neous CCD data corresponding to a sample of 228 elliptical galaxies belonging to 24 clusters of 
galaxies at 0.015 < z < 0.080. We investigate departures from the standard magnitude-isophotal 
size relation, as a function of environmental (cluster-centric distance, local galaxy density) and 
structural (cluster velocity dispersion, Bautz-Morgan type) properties. We find that, for any 
particular galaxy luminosity, the ellipticals in the inner and denser regions of the clusters are 
about 5% smaller than those in the outer regions, which is in good agreement with the finding of 
Strom & Strom (1978d) based on photographic photometry. The null hypothesis (ie., galaxy sizes 
are independent of the cluster-centric distance or density) is rejected at a significance level of 
better than 99.7%. Numericals models of Aguilar & White (1986) predict that tidal stripping can 
lead to changes in the whole structure of ellipticals producing shrinkage and brightening of the 
galaxy, qualitatively consistent with our measurements and also with the findings of Trujillo et al. 
(2002), that more centrally concentrated elliptiials populate denser regions. Our observational 
results can be interpreted as evidence for stripping of stars from ellipticals in the central/denser 
regions of clusters, contributing to the intra-cluster light observed in these structures. 

Subject headings: galaxies: clusters: general- galaxies: clusters: individual: (A3558)- galaxies: 



1. Introduction 

The existence of diffuse optical light in the cen- 
tral region of galaxy clusters is well known (Zwicky 
1951; de Vaucouleurs 1960), but its very low sur- 
face brightness precluded quantitative analysis. 

More recently, due to advances in observational 
and data reduction techniques, accurate measure- 
ments of the total amount of this light and its 
color and spatial distribution became feasible for 
rich clusters (Uson, Bough & Kuhn 1991; Scheick 
k Kuhn 1994; Gonzalez et al. 2000, 2005; Feld- 
meier, CiarduUo & Jacoby 2002; Feldmeier et al. 
2004a). Moreover, single objects, like planetary 
nebulae, red-giant stars or supernovae have been 
detected in the intra-clustcr medium (ICM) (The- 
uns & Warren 1996; Arnaboldi et al. 1996, 2002, 
2003; Feldmeier, CiarduUo & Jacoby 1998; Fergu- 
son et al. 1998; Feldmeier et al. 2004b). 

These studies indicate that between ~10% to 
20% (or higher) of the total stellar light of the 
clusters comes from stars in the ICM. These stars 
have, probably, been stripped from cluster galax- 
ies after galaxy-galaxy (Gallagher & Ostriker 
1972; Richstone 1976) or galaxy-cluster interac- 
tions (Merritt 1984). As a galaxy moves through 
the cluster, it is subject to the gravitational forces 
of their neighbors, as well as of the cluster as 
a whole. As a result of the action of these tidal 
forces, some galaxy stars may be accelerated to ve- 
locities larger than the local escape velocity, being 
removed from the parent galaxy. This process is 
called tidal stripping. In some cases the tidal forces 
may be so strong that the galaxy as a whole is 
disrupted. Indeed, a significant fraction of galaxy 
stars may be dispersed in the ICM after disrup- 
tion of dwarf spheroidals or low surface brightness 
disk galaxies (e.g. Lopez-Cruz et al. 1997), as sug- 
gested by the light plumes observed in clusters 
like Centaurus and Coma (Calcaneo-Roldan et al. 
2000; Korchagin, Tsuchiya & Miyama 2001). 

In the present work we investigate whether stars 

^Visiting Astronomer at the Cerro Tololo Inter- 
American Observatory, National Optical Astronomy Ob- 
servatories, which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc., under a cooperative 
agreement with the National Science Foundation. 

^Thc National Radio Astronomy Observatory is a fa- 
cility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 



stripped of ellipticals can make a significant con- 
tribution to the intra-cluster light (ICL). There 
are some hints that point towards this direction. 
Firstly, ellipticals are basically found in the cen- 
tral, denser, cluster regions, where close encoun- 
ters between galaxies are very common. Their al- 
most radial orbits (Ramirez & de Souza 1998) lead 
them to several close approximations to the cluster 
center, where tidal forces are very strong. The dif- 
fuse light is often distributed like a halo around the 
brightest cluster galaxy (BCG), extending some- 
times beyond ^ 500 kpc. Gonzalez et al. (2005) 
in their study of the surface brightness distribu- 
tion of BCGs in 24 clusters identify a component, 
independent of the BCG, that is possibly due to 
intracluster stars and dynamically linked to the 
cluster as a whole. 

The effects of galaxy-galaxy tidal collisions were 
studied by Aguilar & White (1985, 1986) through 
the use of the impulsive approximation and N- 
body numerical simulations of close collisions be- 
tween galaxies with roughly the same mass and 
with mass distributions consistent with de Vau- 
couleurs and King profiles. They quantified the 
mass and energy changes during such collisions 
and studied the changes of the mass profiles. They 
found that the r^/"* profile is robust, since it is re- 
covered after a collision, but with different param- 
eters. After strong encounters the central surface 
brightness becomes brighter and the effective and 
isophotal radius decrease as mass is stripped out. 
More recently, the high-resolution N-body-|-SPH 
numerical simulations of rich clusters by Willman 
et al. (2004) and Murante et al. (2004) indicate 
that the accumulation of ICL is an ongoing pro- 
cess, linked to infall and stripping events, and that 
the stars composing the ICL are on average older 
than the stars in cluster galaxies. 

Other works have considered the effects of a 
smooth cluster potential on a galaxy. Using the 
impulsive approximation (Merritt 1984; Mamon 
2000), or high resolution numerical simulations 
(Ghigna et al. 1998), it can be shown that clus- 
ter galaxies are truncated at the tidal radius: 

rud = C (^^) rp, (1) 

where Ug and ad are the velocity dispersions of 
the galaxy and the cluster, respectively, rp is 
the radius of maximum approximation between a 
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galaxy and the cluster center, and C is a non- 
dimensional constant that depends of several fac- 
tors, like galaxy orbits, and is close to unity. This 
process is often called tidal truncation. 

It is not clear yet what process is more effec- 
tive in removing stars from cluster galaxies. Tidal 
truncation is indeed more efficient in removing 
halo stars, since rud is larger than the optical ra- 
dius for the majority of the galaxies. Anyway, 
as pointed out by Gnedin (1999), mass is not 
only lost by instantaneous stripping after cluster- 
galaxy shocks, but also by secular tidal heating. 

The first observational detection of tidal strip- 
ping of stars from elliptical galaxies came from 
Strom & Strom (1978a, b,c,d), with a sample of 400 
ellipticals and SO's located in 6 different clusters, 
deeply imaged with photographic plates. These 
authors found that ellipticals in the center of rich 
clusters have isophotal radii 10-30% smaller at a 
given magnitude than ellipticals in poor clusters 
or in the outskirts of rich ones, estimating that up 
to ^ 37% of the luminous mass of the stripped 
ellipticals have been lost to the ICM. 

Despite of the great impact of this work, it 
has been severely criticized, mainly because it was 
based on photographic photometry (e.g. Dressier 
1984). Giuricin et al. (1989) questioned the sta- 
tistical analysis done by the Strom & Strom, 
and using photo-electric and CCD data of ~ 160 
elliptical galaxies from the sample of Burstein 
et al. (1987), did not find any dependence of 
the magnitude-size relation properties of ellipti- 
cal galaxies with the environment. In this paper 
we will revisit this topic, using CCD data homo- 
geneously collected and reduced. 

The layout of the paper is as follows. In Section 
2 we describe the observational material and the 
cluster sample, we present our method for the ex- 
traction of galaxy photometric parameters, as well 
as the selection criteria of the final galaxy sam- 
ple. In Section 3 we present our results on galaxy 
sizes, that are discussed in Section 4. In Section 5 
we summarize the main conclusions of this work. 
Throughout this paper, we adopt, when necessary, 
a ACDM cosmology with Hq = 70 hjo km s~^ 
Mpc-\ = 0.3 and = 0.7. 



2. The data 

2.1. The cluster sample and the observa- 
tions 

The galaxy sample discussed in this work was 
extracted from the set originally obtained as part 
of the Ph.D. thesis of Daniel A. Dale on pecu- 
liar motions of clusters with z < 0.1 (Dale et al. 
1997, 1998, 1999a,b; Dale et al. 1999c). The 
original cluster sample was selected from nearby 
{cz < 18000 km s"^) objects of the Abell, Cor- 
win & Olowin (1989) cluster catalogue, in order 
to cover as much and uniformly the whole sky 
area as possible. Although the cluster selection 
criteria and imaging characteristics were chosen 
to optimize the peculiar motion study, the result- 
ing database is useful for many other studies, since 
it was homogeneously observed and reduced, has 
high photometric accuracy, and the cluster sample 
is representative of the low redshift cluster popula- 
tion (Cypriano et al. 2001). For the present study, 
only the Southern part of this sample has been 
used. The clusters, and some of their properties 
that are relevant for this work, are presented in 
Table 1. 

The observational material consists of several 
Kron-Cousins I band images obtained with the 
0.9m CTIO telescope. The details of the obser- 
vations are discussed elsewhere (Dale et al. 1997, 
1998) and here we only summarize them. The de- 
tector used was the 2kx2k Tek2k No. 3 CCD, with 
a scale of 0.4 arcsec per pixel, resulting in a field 
of 13.5' X 13.5' per image. The exposure time was 
600 seconds in all cases. The images reach 23.0 to 
23.8 I mag arcsec"^ at the 1 a level over the back- 
ground, with a median of 23.6 1 mag arcsec"^. The 
median seeing of the images is 1.52" and ranges 
from 1.16" to 2.80". 

These images, consisting of several pointings 
per cluster, in general do not uniformly cover a 
cluster, since they were taken in regions near spiral 
galaxies, although the central region of the clusters 
is always covered. For each cluster, the fraction of 
the projected area inside the R200 radius (see sec- 
tion 2.4) covered by the imaging is also shown in 
Table 1. The actual distribution of pointings is 
presented in Dale et al. (1997, 1998, 1999b) 
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Table 1 
Cluster sample 



Name 


z 


cl 


B-M Type 


Totcil cirecL 


Fraction of area 


N 






(km s"^) 




(squcire degree) 


(%) 






(2) 


f3) 


f4) 


(5) 


(6) 


(7) 


A85 


0.0555 


941 


I 


0.05 


4.8 


12 


A114 


0.0582 


911 


II 


0.39 


44.8 


15 


A119 


0.0449 


685 


IITII 


0.05 


5.3 


6 


A496 


0.0326 


686 


I 


0.28 


16.6 


19 


A2670 


0.0766 


873 


HI 


0.05 


11.4 


7 


A2806 


0.0265 


433 


HI 


0.26 


34.0 


4 


A2877 


0.0251 


1026 


I 


0.45 


10.2 


7 


A2911 


0.0800 


546 


HI 


0.32 


12.8 


2 


A3193 


0.0346 


638 


I 


0.23 


21.3 


5 


A3266 


0.0586 


1131 


I-II 


0.18 


15.0 


17 


A3376 


0.0456 


641 


I 


0.10 


14.2 


8 


A3381 


0.0372 


414 


I 


0.09 


39.3 


2 


A3389 


0.0259 


511 


II-III 


0.10 


5.8 


1 


A3395 


0.0496 


823 


II 


0.10 


9.9 


10 


A3407 


0.0415 


504 


I 


0.13 


29.3 


3 


A3408 


0.0405 


900 


I-II 


0.20 


12.4 


4 


A3528 


0.0538 


955 


II 


0.31 


26.7 


16 


A3558 


0.0473 


935 


I 


0.56 


41.4 


35 


A3571 


0.0384 


969 


I 


0.05 


2.1 


6 


A3574 


0.0152 


447 


I 


0.09 


2.2 


1 


A3656 


0.0195 


366 


I-II 


0.19 


15.5 


5 


A3667 


0.0549 


987 


I-II 


0.44 


33.9 


19 


A3716 


0.0454 


842 


I-II 


0.40 


33.4 


13 


A3744 


0.0386 


624 


II-III 


0.31 


37.9 


7 


A4038 


0.0302 


826 


III 


0.33 


11.3 


4 



^Wc got rcdshifts and velocity dispertions form Abell cluster rcdshift compilation of 
Andcrnach & Tago (2005, in prep.), see Andernach, Tago (1998) for a description. 



Note. — (1) name of the cluster in the Abeh catalog (Abeh, Corwin & Olowin 1989); (2) 
redshift; (3) radial velocity dispersion; (4) Bautz-Morgan Type (Bautz & Morgan 1970); 
(5) total imaged area; (6) The ratio of the area covered by the imaging over the whole 
area inside the i?200 radius ; (7) number of ellipticals in the sample. 
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2.2. Determination of photometric param- 
eters 

Most of the images used here were obtained 
under good photometric conditions. The abso- 
lute cahbration of the magnitudes was done by us- 
ing standard photometric stars of Landolt (1992). 
The photometric zero-point cahbration could be 
determined with a median accuracy of 0.018 mag, 
never larger than 0.031 mag. 

In this work we have used an isophotal radius 
as an estimator of the galaxy size because its log- 
arithm is very well correlated with the magni- 
tudes (Spearman rank correlation coefficient = 
—0.99). On the other hand, the log of the radius 
that encloses half of a galaxy light (the effective 
radius r^) shows a significantly poorer correlation 
with magnitudes (r^ = —0.86). The Petrosian ra- 
dius (Petrosian 1976) is also less well-correlated 
with magnitudes (r^ — —0.84), which is is not un- 
expected, since elliptical galaxies spanning a range 
of 5 absolute magnitudes have different profiles 
(Caon, Capaccioli & D'Onofrio 1993). 

The isophotal quantities are defined at the same 
isophotal limit in the galaxy rest frame (fiHm)- 
The corresponding isophotal threshold (/Xt), in the 
observer rest frame, is given by the following ex- 
pression: 

+ Ai + 10\og{l + z) + k{z) (2) 

where A/ is the Galactic absorption in the I band, 
101og(l -I- z) is the correction due to the cosmo- 
logical dimming, and k{z) is the /c-correction. 

The values for Galactic absorption adopted here 
are those given by Schlegel, Finkbciner & Davis 
(1998), and range from 0.032 to 0.363 for the clus- 
ters in the sample. The fc-corrections are from 
Poggianti (1997), assuming her model for ellipti- 
cal galaxies, and range from 0.006 for the near- 
est cluster (A3574, z=0.014) to 0.035 mag for the 
farthest one (A2670, z=0.076). The cosmological 
dimming factor ranges from 0.140 to 0.318 mag. 
The median value of fit — Mzim is 0.32 mag. 

The surface photometry was performed with 
the task elUpse, of IRAF/STSDAS ^ The output 
of this program is a list of several parameters, such 

^IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of Uni- 



as ellipticities, local intensities, magnitudes, etc., 
as a function of the semi-major axis value. Using 
as a radius the geometrical average of the semi- 
major and semi-minor axis, the so-called equiva- 
lent radius, the isophotal radius (r^so) is estimated 
from the galaxy isophotal brightness profile as the 
radius corresponding to fit , and then the isophotal 
magnitude {rriiso) inside Viso is determined. Using 
the equivalent radius instead of the semi-major 
axis we remove possible bias due to different el- 
lipticities between galaxies. 

A de Vaucouleurs (1948) and an exponential 
profile were also fitted to the isophotal bright- 
ness profiles. Following the procedure described 
in Scodcggio, Giovanelli & Haynes (1998), we re- 
moved from the fitting an inner region, with radius 
equal to two times the FWHM of stellar sources. 
These fittings were used as part of the selection of 
elliptical galaxies (see Section 3.3). 

The median of the formal errors are 0.02 in 
magnitudes and 3.3% in r^so, where errors in riso 
have been calculated as: 

dv 

o-(riso) = o-A"^- (3) 

with (T^ being the error of the corresponding sur- 
face brightness. 

A check of the internal consistency of these er- 
rors can be obtained from a sample of 16 galax- 
ies that are present in two different images. The 
rms differences between isophotal magnitudes and 
isophotal radii measured in different images are 
0.02 and 2.7%, respectively, both in agreement 
with the formal errors. 

2.3. Sky subtraction 

Since the sky subtraction is critical in the kind 
of study presented here, the procedure of sky de- 
termination is now explained in detail. 

Firstly, the software SExtractor (Bertin & 
Arnouts 1996) is used to determine the average 
sky value for the entire field. Then the galaxy 
image is modeled using the IRAF task ellipse and 
removed from a sky subtracted image. After, we 
estimate the residuals in the local sky, by examin- 
ing the residuals in a square region centered in the 

versities for Research in Astronomy, Inc. (AURA) under 
cooperative agreement with the National Science Founda- 
tion. 
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subtracted galaxy, with a side size equal to at least 
two times the galaxy major-axis at fxt- The size 
of this region was determined after several tests 
and proved to be convenient for the estimation 
of the local background. The new local sky level 
is obtained by removing the mean value of the 
residuals from the previously adopted background 
level. 

After, the procedure of surface photometry and 
image subtraction is done again, using this new 
value for the value of the local sky. The whole pro- 
cess is repeated until no significant local sky vari- 
ations are found. Generally, two iterations were 
enough. 

This process starts with the brightest galaxy 
in the field and continues in order of increasing 
magnitude, to avoid luminosity contamination of 
the galaxy by the envelope of the BCG and other 
bright galaxies in the cluster field. 

2.4. Selection of the galaxy sample 

The selection of elliptical galaxies for the anal- 
ysis was done following several steps. 

A first selection was done based on the galaxy 
morphology and the form of its light profile. Ini- 
tially, all galaxies that are obvious spirals were 
removed from the sample. Then, by comparing 
galaxy light profiles with exponential and de Vau- 
couleurs (1948) profiles, it was possible to remove 
from the sample of elliptical galaxies some appar- 
ently featureless disk-like galaxies. 

After, the morphological classification was 
checked with previous classifications using NED 
^. Nearly half of the galaxies in the sample had 
previous classifications. About 30% of them are 
classified as SO and were removed from the sam- 
ple. Only a few cases of early-type spirals classified 
as ellipticals were found. Galaxies with redshift 
inconsistent with those of the clusters were also 
removed. 

Due to the presence of a disk component, SO 
galaxies tend to have larger sizes, at a given mag- 
nitude, than ellipticals. However, since the radial 
distribution of SOs and ellipticals in clusters does 



^This research has made use of the NASA/IPAC Extragalac- 
tic Database (NED) which is operated by the Jet Propul- 
sion Laboratory, CaUfornia Institute of Technology, under 
contract with the National Aeronautics and Space Admin- 
istration. 



not strictly match each other, there is the possi- 
bility that misclassification of lenticulars in ellip- 
ticals might produce a radial trend that could bias 
our results, because the radial distribution of the 
former is slightly more extended than the latter 
(Whitmore, Gilmore & Jones 1993). This is diffi- 
cult to evaluate because we do not know, a priori, 
neither the number of misclassifications nor their 
radial distribution. However, we expect that they 
would be more common in the inner than in the 
outer parts of the clusters, due to the difficulty 
introduced by the cluster diffuse light in the de- 
tection of faint disks. If this is indeed the case, 
the bias introduced would only dilute the radial 
trends discussed here. 

To ensure a high photometric quality of the 
data, only galaxies with isophotal radius larger 
than two times the FWHM of point sources plus 
0.8" (2 pixels) were included in the sample. Addi- 
tionally, images where the chosen isophotal limit is 
within 1(7 of the sky brightness were excluded from 
the analysis. The faintest galaxies of the sample 
have 17.5 I mag. 

Objects with signs of contamination by a 
nearby object, which could be easily recognized 
on the light profiles and in the residuals that re- 
main after the subtraction of the galaxy, were also 
removed from the sample. Objects with strong 
gradients in the local background, precluding an 
accurate sky subtraction, were also removed. 

As a final criterion, only galaxies with cluster- 
centric distances lower than R200 were included in 
the analysis. This is the radius where the mass 
density is ~200 times larger than the critical den- 
sity and is a good approximation to the virial ra- 
dius. This radius R200 is estimated using the ex- 
pression 



.R200 — 



10H{z) 



(4) 



(Carlberg et al. 1996) where ad is the one- 
dimensional velocity dispersion of the cluster and 
H(z), that is a function of the cosmological param- 
eters, is the Hubble factor at the cluster redshift. 

The size of the sample depends of the isophotal 
limit chosen. Using a fiducial isophotal limit in the 
galaxy rest frame of = 22.75 I mag arcsec"^ 
, which avoids significant incompleteness, the full 
sample contains 228 elliptical galaxies, 172 (75%) 
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of them with spectroscopic confirmation of cluster 
membership. The number of eUipticals per cluster 
in the final sample is presented in Table 1. 

At this point it is important to warn the reader 
that the galaxy sample selected for this work is not 
complete in any sense, due to the way the fields 
of the clusters were imaged and also because we 
have rejected a number of objects for which we 
were not able to perform the surface photometry 
properly, mainly because of light contamination 
by a neighbor. This problem is more severe in 
the central parts of the clusters, since there the 
surface density of objects is large. Both factors of 
incompleteness can change the number of objects 
as a function of the radius, but it is not expected 
that they can lead to any bias that can explain the 
trends identified in this work (see next section), 
because they are completely independent of the 
photometric parameters of the galaxies. 
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3. Results 

3.1. The magnitude-size relation for clus- 
ter ellipticals 

In Figure 1 we plot the observed relationship 
between the absolute isophotal magnitude and the 
isophotal radius (at the fiducial isophotal limit 
in the galaxy rest frame of fJ-Hm = 22.75 I mag 
arcsec"^), as well as standard least-square fits of 
linear and quadratic models. 

It can be appreciated in this plot that the linear 
model (dashed line) fits well the data, but with 
some systematic residuals in both the faint and 
bright extremes. However, these residual are min- 
imized by the use of a quadratic model (continu- 
ous line). Actually, a likelihood ratio test indicates 
that a quadratic model is significantly more reli- 
able than a linear model. Higher order models do 
not improve statistically the likelihood of the fit. 

Hereafter, <riso>{M) will be used to denote 
the isophotal radius of a galaxy that, given its I 
magnitude, obeys the magnitude-size relation rep- 
resented by the continuos line in Figure 1. We 
define 



where ri^o is the actual isophotal radius, as a 
convenient way of measuring the deviation of the 
isophotal radius of a galaxy from the value that is 
expected given its luminosity. Note that rj < 1 rep- 



Fig. 1. — Size-Magnitude relation in the I band 
for cluster ellipticals. The isophotal limit where 
radii and magnitudes are measured, in the galaxy 
rest frame, is num = 22.75 1 mag arcsec"^. The 
straight dashed line and thr continuous line are 
least-squares fit to the data using a linear and a 
quadratic models respectively. 

resents an effective "shrinkage" of a galaxy. The 
mean uncertainty in rj is 0.037, taking into account 
the error in the measured isophotal radius and in 
the fitted radius-magnitude relation. 

We present in Figure 2 the distribution of rj for 
all galaxies in our sample. Its average value is 
1.005 with a standard deviation of 0.10 (or 10%). 

The observational errors in r.iso can account for 
one third of this scatter and the errors on the mag- 
nitudes are negligible (formal errors in Mj are typ- 
ically 0.02 mag). Indeed, most of the variance in 
77 is due to cosmic scatter in the magnitude-size 
relation. 

We verified that the image quality is not affect- 
ing rj measurements. The Spearman rank correla- 
tion coefficient between rj and the seeing (FWHM) 
of the images is 0.1, indicating an inexistent or 
very low dependence between these two quanti- 
ties. 
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Fig. 2. — Distribution of values of the ratio be- 
tween the measured isophotal radius of a galaxy 
and that expected for a galaxy with the same mag- 
nitude, T], for the galaxies of the sample. 



3.2. Galaxy sizes as a function of cluster- 
centric distances 

In the seventies, Strom & Strom claimed that 
they had found evidence for tidal stripping in 
galaxy clusters based on an analysis of the isopho- 
tal radius of cluster galaxies and a function of 
its cluster-centric radial position (Strom & Strom 
1978d). We carry out a similar analysis in this 
section using the cluster galaxies of our sample. 

Since the number of data points per cluster is 
generally low, we have decided to use an ensemble 
technique, where the data of all clusters are ana- 
lyzed together. For that, the cluster-centric dis- 
tances of the galaxies, R, were normalized by the 
value of the i?200 of the cluster. Figure 3 presents 
the ratio as a function of the normalized cluster- 
centric distance R/R200 for our sample. 

To investigate radial trends in the optical ex- 
tent of elliptical galaxies, we examine whether rj 
varies with i?/i?2oo- Indeed, by dividing the sam- 
ple with respect to the median of R/R200 (0.202), 
we find that galaxies in the group near the center 
((i?/i?20o) 0.11) have 77 = 0.980±0.009, whereas 
those in the outer group ((i?/i?20o) = 0.43) have 



Fig. 3. — Mean value of 77 as a function of the 
normalized projected radial distance i?/i?2oo- The 
main plot shows bins with approximately the same 
number of data points. The horizontal error bars 
are the 1-sigma dispersion of the data and the ver- 
tical ones the error of the mean (a/^/N). The 
small panel shows the individual data points. On 
the upper right corner of the small panel is shown 
an average rj error bar. 

f] ~ 1.030 ± 0.009. The mean difference in 77 be- 
tween the two groups is (5. Oil. 3)%, ruling out at 
a 3.8cr (99.99%) confidence level (C.L.) the null hy- 
pothesis that sizes are independent of the cluster- 
centric distance. 

Comparing the medians, which are more robust 
estimators of location than the means, galaxies of 
a given luminosity in the outer group are 5.1% 
larger than those in the inner regions of the clus- 
ters. 

It is interesting to verify how this result changes 
as a function of fium ■ In Table 2 we show the dif- 
ference in 77 (for medians and means as estimators) 
of galaxies in these two groups of projected radial 
distances, for different isophotal limits. For both 
estimators, this difference tends to increase toward 
fainter values of fium , although this gradient is well 
inside the errors. 

The results of a Student's t-tcst for the statisti- 
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cal significance of the difference of the mean value 
of r] for the two radial groups arc also shown in 
Table 2 and indicates that this difference is signif- 
icant for all values of Hum ■ Note that the difference 
between the medians is almost always larger than 
that between the means. 

The trend suggested by the data in Table 2 is 
that the differences between rj of the radial groups 
increases for fainter isophotal limits. This is par- 
ticularly clear for the difference between medians. 
A word of caution here is that the significance of 
this gradient is low due to the errors in these dif- 
ferences. We show in Figure 3 mean values of 77 
for four radial bins with approximately the same 
number of data points in each one. The figure sug- 
gests that 7] increases monotonically with radius, 
possibly reaching a plateau at R/R200 ^ 0.4. 

Unfortunately, most of the clusters do not 
have a sufficiently large number of galaxies in 
our sample to allow separate analyses. How- 
ever, Abell 3558 is an exception. This cluster, 
in the center of the Shapley super-cluster, is very 
rich and has 36 galaxies in our sample, allowing 
a statistical analysis. The 77-radius relation for 
Abell 3558 is very significant, as shown in Figure 
4. The difference between mean values of 77 be- 
low and above R/R200 = 0.202 (0.68 h^^ Mpc) is 
(9.8 ±3.6)%. 

In Table 3 we show the difference of rj between 
the inner and the outer radial groups of galaxies 
for the 5 clusters which have at least 15 galax- 
ies in our sample. These clusters are among the 
richest in Table 1 . In four out of the five clusters 
the galaxies in the outer group tend to be larger 
than those in the inner group, although for one of 
them (A3266), the difference is not really signifi- 
cant. For A3528 the difference is in the opposite 
direction but, again, the significance is low. 

3.3. Galaxy sizes as a function of the local 
galaxy density 

We are also interested on verifying whether 
r] varies with the local galaxy density. Since 
our imaging covers the clusters only partially, to 
avoid incompleteness in the spatial sampling we 
have used Super COSMOS data (Hambly, Irwin 
& MacGillivray et al. 2001) to obtain a catalog of 
the galaxies projected on to the cluster fields. We 
have selected only galaxies brighter than R = 16.5. 



A3558 



1.10 - 



1.05 



1.00 



0.95 - 



0.90 




0.2 0.4 0.6 0. 



Fig. 4. — Distribution 77 as a function of the 
normalized projected radial distance i?/i?200 for 
galaxies in Abell 3558. 



The local surface density associated to each 
galaxy in our sample has been computed from the 
projected distance to its sixth closest Super COS- 
MOS neighbor (Re), using the estimator of Caser- 
tano & Hut (1985) 



(6) 



nRl 



This quantity should be corrected by contamina- 
tion of foreground and background galaxies. 

This field contamination was estimated as the 
median surface density of galaxies in a ring at 
4.5±0.2 R200 of the cluster center. However, this 
correction is quite sensitive to the variance of the 
field galaxy counts, often produced by large scale 
structures near the cluster or in its line of sight. 
The cluster A3558, in the center of the Shapley 
super-cluster, illustrates well this problem. De- 
spite its richness (Abell richness class 4), the den- 
sities estimated here are among the lowest of the 
whole sample because the background counts are 
strongly enhanced due to the presence of several 
clusters in the region. For this reason, we also 
analyzed the data without any background sub- 
traction. 
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Table 2 

Variation of differences of rj for galaxies in the outer and inner radial groups, as a 

FUNCTION OF THE LIMITING ISOPHOTE 





Number 


Difference between 


Student's t-test 


(mag arcsec"^) 


of galaxies 


medians (%) 


means (%) 


t 


Prob. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


21.50 


151 


5.4 


4.7 ± 1.1 


4.4 


100.00% 


21.75 


167 


4.9 


5.0 ± 1.1 


4.5 


100.00% 


22.00 


179 


6.1 


4.8 ± 1.2 


4.2 


100.00% 


22.25 


191 


5.4 


4.6 ± 1.3 


3.7 


99.97% 


22.50 


210 


4.8 


5.0 ± 1.2 


4.0 


99.99% 


22.75 


228 


5.1 


5.0 ± 1.3 


3.8 


99.98% 


23.00 


200 


8.2 


7.0 ± 1.6 


4.5 


100.00% 



Note. — (1) Isophotal limit in the I-band in the galaxy rest frame; (2) Number 
of galaxies; (3) Differences between the medians; (4) Differences between the 
means and their errors; (5) Value of t from the Student's <-test; (6) Probability 
that the observed difference between the two radial groups will not be caused by 
statistical fluctuations. 



Table 3 

Variation of differences of ry for galaxies in the outer and inner radial groups for 

INDIVIDUAL clusters 



Cluster 


Number 


Difference between 


Student's t-test 




of galaxies 


medians (%) 


means (%) 


t 


Prob. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


A3558 


35 


13.6 


9.8 ± 3.6 


2.5 


98.3% 


A3667 


19 


11.1 


10.2 ± 3.8 


3.0 


99.3% 


A496 


19 


14.4 


15.7 ± 4.6 


3.5 


99.7% 


A3266 


17 


1.9 


1.6 ± 2.9 


0.5 


36.0% 


A3528 


16 


-9.7 


-6.9 ± 6.3 


-1.2 


-75.2% 



Note. — Only clusters with more than 15 galaxies on the sample 
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Figure 5 shows the relation between rj and the 
logarithm of the corrected projected number den- 
sity. The figure suggests that there is a tendency 
for rj to decrease with increasing surface density. 

The difference of 77 between galaxies in envi- 
ronments with densities larger and smaller than 
the median for the background subtracted sample 
is 3.1% (median) and 3.3 ± 1.3% (mean; 98.6% 
C.L.). The values for the non-background sub- 
tracted densities are similar: 3.3% (median) and 
4.0 ± 1.3% (mean; 99.7% C.L.). The median value 
of E for the background subtracted sample is 51.3 
Mpc"2 and 84.3 Mpc"^ foj. ^i^^^ 

non-background 

subtracted one. 



1.04 



1.02 



1.00 



0.98 



0.96 



0.94 



log E (gal. Mpa-') 



1 2 

log S (gal. Mpc"^) 



Fig. 5. — Distribution of 77 as a function of the log- 
arithm of the projected number density of galax- 
ies, E. The solid circles shows averages of bins 
with similar number of data points for the back- 
ground subtracted values of the density. The 
open circles correspond to the non background- 
subtracted densities. The horizontal error bars are 
the 1-sigma dispersion of the densities in each bin 
and the vertical ones the error of the means. In 
this plot negative densities are represented as the 
logarithm of the absolute density value times -1. 



3.4. Galaxy sizes as a function of cluster 
properties 

In order to understand the origin of the strip- 
ping mechanism it is important to verify whether 
and how the isophotal diameters of the galaxies 
depend on cluster properties. Here we will ana- 
lyze two of them: velocity dispersion and Bautz- 
Morgan type that, in principle, are approximately 
related to mass and evolutionary phase, respec- 
tively. In Figures 6 and 7 we show the dependence 
of 77 as a function of these parameters. 

1.05 



1.00 



0.95 



400 600 aoo 1000 I 

Vel, dispersion (km a"') 




I & 



400 



600 



800 



1000 



Vel. dispersion (km s ') 



1200 



Fig. 6. — Distribution of 77 as a function of the 
dispersion of velocity of the cluster. The circles 
represents the mean value of 77 with its error of 
the mean; the horizontal bars represent the inter- 
val of ad associated with each bin. Filled circles 
represent mean values for all galaxies in the bin, 
whereas the open ones contain only galaxies with 
cluster-centric distances less than 0.202 R2oo- For 
clarity, we shifted horizontally the open circles by 
20 km s-i. 

At a first glance neither of these figures dis- 
plays a clear correlation, as those found for the 
cluster-centric distance and the local density. In 
Figure 6 it can be seen that galaxies in clusters in 
the smallest velocity dispersion bin, present values 
of 77 significantly smaller than those in the other 
clusters. 
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1.05 



1.00 




I-[[ II II-[II 

Baulz-Morgan Type 



0.95 



_L 



I-II II II-III 
Bautz — Morgan Type 



III 




Fig. 8. — {Left panel) Mean value of the magni- 
tude as a function of the cluster-centric distance. 
{Right panel) Mean value of 7y as a function of the 
magnitude. In both panels, the bins represent ap- 
proximately the same number of data points. The 
horizontal error bars are the 1-sigma dispersion 
of the data and the vertical ones the error of the 
mean {a/VN). 



Fig. 7. — Same as Figure 6, but as a function of 
the Bautz-Morgan Type. 

Figure 7 shows how 77 varies with the morpho- 
logical type of the clusters. Again we cannot see 
any trend here. What is noteworthy in this figure, 
however, is that elliptical galaxies (excluding the 
BCG) in the inner region of cD clusters (Bautz- 
Morgan type I) tend to be smaller than those in 
other cluster types, especially those closest to the 
cluster center. Actually, the inner galaxies of B- 
M type I clusters are the smallest in the present 
sample. A similar result was obtained by Sandage 
& Hardy (1973), who found that elliptical and 
lenticular galaxies in type I clusters are dimmer, 
at given radius, than in type II clusters. 

3.5. Is there a dependence of rj with mag- 
nitude? 

An important check is to verify whether the 
effects found here are an artifact of the fitted 
magnitude-size relation or the result of luminosity 
segregation. We address this issue in Figure 8. 

The left panel of Figure 8, presenting the lu- 
minosity distribution as a function of the cluster- 
centric radius, shows a hint that the brighter 
galaxies tend to be found in the outer regions. 
The difference considering only two radial bins, 
however, is rather insignificant, 0.27 ± 0.18 mag. 



Moreover, the right panel reveals that there is not 
a monotonic dependence between magnitude and 
77. For instance, there is no significant relative size 
difference between the bright and faint half of the 
sample, 0.8±1.4%. In fact, this plot shows a dif- 
ference between rj values at the extremes and the 
center of the magnitude distribution. An ispec- 
tion of Figure 1 shows that this is produced by the 
residuals of the fit of the magnitude-size relation 
with our quadratic model. Rather than indicating 
segregation (or anti-), the trends in Figure 8 are 
probably reflecting features of the galaxy sample, 
which is not complete in any sense (c.f. Sect. 2). 

For a double check, we repeated our analy- 
sis considering only galaxies with magnitudes be- 
tween -23.5 and -20.5. Within this interval, the 
linear and quadratic fits to the magnitude-size re- 
lation are nearly identical. We obtained, by do- 
ing that, essentially the same results, although 
they are noiser due to the smaller number of data 
points. We thus conclude that all the results found 
in this work, regarding the relative size of ellip- 
tical galaxies in clusters, are not artifacts of the 
magnitude-size relation fit or a function of the 
galaxy luminosity. 
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4. Discussion 

4.1. Comparison with previous work. 

Our main result, the shrinkage of galaxies in- 
habiting the inner parts of the clusters, is in qual- 
itative agreement with the results of Strom & 
Strom (1978d). These authors found shrinkage 
factors of 10-30% in images reaching isophotal lev- 
els as faint as 26 R mag arcsec"^, or ^^25. 3 I 
mag arcsec"^, assuming (R-I) = 0.7, which is typ- 
ical of ellipticals at z~0 (Fukugita, Shimasaku & 
Ichikawa 1995). This is nearly 2.5 mag arcsec"^ 
fainter than our fiducial isophotal limit, 22.75 I 
mag arcsec"^, precluding a quantitative compar- 
ison between our results in Table 2 with those of 
Strom & Strom (1978d), because the reduction of 
the isophotal radius may depend on the isophotal 
level. 

Giuricin et al. (1989) have not found any en- 
vironmental dependency on the sizes of ellipti- 
cal galaxies in their sample. Their photometry 
reaches 25.0 B mag arcsec"^, which is roughly as 
deep as our fiducial value for fium (assuming (B- 
I) « 2.27 for ellipticals; Fukugita, Shimasaku & 
Ichikawa 1995). However, they warned in their 
conclusions that a shrinkage of the isophotal radii 
of a few percent is consistent with their results, 
given their observational uncertainties. Actually, 
Giuricin et al. (1989) report that their typical rms 
errors in logr^so are around 0.04 (or 9.6%), and, in 
the best cases (CCD data) they are between 0.02 
to 0.03 (4.7% to 7.2%), weU above the rms er- 
ror of the log of the isophotal radius in this work, 
which is 0.015 (3.5%), as explained in Section 2. 
Additionally, the present work uses the ensemble 
technique, that allows putting together data from 
different clusters, while Giuricin et al. (1989) di- 
vided their sample in sub-samples with small num- 
ber of data points, leading to results with smaller 
statistical significance. 

Trujillo et al. (2002) have found a related re- 
sult through the analysis of light concentration for 
the Virgo and Coma clusters, indicating that more 
centrally concentrated elliptical galaxies tend to 
inhabit the denser regions of galaxy clusters, which 
is in agreement with the results presented here. 
When interpreting their results, these authors sug- 
gest that the origin of such effect would be galaxy 
mergers, rather than tidal encounters because, 
quoting Strom & Strom (1978d), the latter would 



affect only the stars in the outer halo, leaving the 
cores untouched. However, the numerical simu- 
lations of Aguilar & White (1986) convincingly 
demonstrate (c.f. Figure 2 of this paper) that 
tidal shocks do lead to changes in the whole galaxy 
structure, not just in the outer parts: "...strong 
collisions produce a shrinkage in re and a bright- 
ening of fie, whereas weak collisions have the op- 
posite effect" . Here le and fXe are the parameters 
of the de Vaucouleurs profile: the effective radius 
and the surface brightness inside Ve- 

Although tidal stripping may not be the only 
mechanism originating the observed shrinkage of 
ellipticals, this same mechanism has been vali- 
dated through numerical simulations (Willman et 
al. 2004; Richstone 1976) to explain the observed 
intracluster light in clusters of galaxies. One of the 
findings of our work, that, for cD clusters (B-M 
type I), the shrinkage of the central ellipticals is 
more pronounced than in the other cluster types, 
is consistent with the early Richstone (1976) claim 
that cD halos can be formed by stars stripped out 
from cluster galaxies by tidal shocks, also in agree- 
ment with the very recent work of (Gonzalez et al. 
2005). 

4.2. Tidal stripping mechanism 

Finally, we discuss how can it be decided ob- 
servationally whether tidal truncation (cluster- 
galaxy collisions) or coUisional stripping (galaxy- 
galaxy encounters) is the dominant mechanism for 
the shrinkage of the ellipticals. A way to distin- 
guish between these alternatives , is by examining 
how the stripping efficiency (and thus the shrink- 
ing of a galactic radius) depends on the cluster 
velocity dispersion. In the first case, as shown in 
equation 1, the efficiency should be proportional 
to ad, because the tidal radius depends on c^^; 
while, in the second case, the dependency should 
be the opposite, since lower velocity dispersions 
leaves room for longer lived, and hence stronger, 
galaxy-galaxy shocks. Such test applied to our 
present data (Section 3.4), unfortunately does not 
show a clear correlation between galactic sizes and 
ad- An analysis of a much larger sample, and 
preferably deeper, is necessary to determine which 
is the dominant mechanism producing the tidal 
stripping of ellipticals in clusters of galaxies. 
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5. Conclusions 

We have analyzed a sample of 228 elliptical 
galaxies belonging to 24 clusters to look for ev- 
idence of tidal stripping and to put constraints 
on the mechanisms behind star removal in cluster 
galaxies. The stripping was quantified by examin- 
ing departures of the magnitudc-isophotal radius 
relation for cluster ellipticals, measured by the pa- 
rameter rj- the ratio between the measured isopho- 
tal radius of a galaxy (at 22.75 I mag arcscc"^ in 
the galaxy rest frame) and the isophotal radius 
expected for a galaxy with the same magnitude. 

Our main conclusion is that cluster ellipticals of 
a given magnitude in inner/denser cluster regions 
are smaller than those in the outer regions by a 
factor of the order of 5%. The shrinkage factor 
tend to decrease with increasing distance from the 
center of the cluster and with decreasing galaxy 
number density. Although the amount of galaxy 
stellar mass lost to the ICM remains unknown, it is 
probably of the order of a few percent of the parent 
galaxy mass within the optical radius. Thus, it is 
highly probable that the stars lost by the central 
cluster ellipticals are a significant contributor to 
the diffuse light observed in clusters of galaxies. 
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